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Abstract 

In a recent paper [J. Stat. Mech. P07006 (2004)], E.G.D. Cohen and David Mauzerall (CM) 
have argued that the derivation of the nonequilibrium work relation given in [C. Jarzynski, Phys. 
Rev. Lett. 78, 2690 (1997)] is flawed. Here I attempt to answer their criticisms, both by presenting 
a detailed version of that derivation and by addressing specific objections raised by CM. The 
derivation presented here is in fact somewhat stronger than the one I gave in 1997, as it does not 
rely on the assumption of a weak coupling term between the system of interest and its thermal 
environment. 
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In a recent paper Cohen and Mauzerall (CM) have raised questions about the va- 
lidity of the nonequilibrium work relation, EqQ below, which relates the external work W 
performed on a system during a nonequilibrium process, to the free energy difference AF 
between two equilibrium states of the system. This prediction has been derived by various 
means in e.g. Refs. l1 lllA and has been confirmed in an experiment 

performed by Liphardt et alVl , along lines suggested by Hummer and Szaboj^, involving 

ri 

the forced unfolding and refolding of a single strand of RNA. [13| Two recent papers review 
theoretical, computational, and experimental aspects of the nonequilibrium work relation 

are aimed at a derivation 



of EqJTJthat I published in 1997 [2j, CM also suggest that other derivations (in particular 
that given by Crooks in Ref. jj]) are relevant only for near-equilibrium processes, and that 
the experiment of Ref. ^ cannot be viewed as a confirmation of Eq^ 

Some of the issues raised by CM pertain to aspects of Eq^that are counter-intuitive, or 
at least sufficiently unusual to arouse justifiable skepticism. For instance, while the right 
side of Eq^is familiar enough - it is a ratio of partition functions - the left side is not. The 
quantity inside the angular brackets is constructed by combining two values - the work W 
performed on a system that is driven out of equilibrium, and the temperature T of the initial 
equilibrium state of the system - into something that looks like a Boltzmann factor, namely 
e _/3W/ (where 1//3 — kgT). To CM this is an ad hoc, unjustified construction, and I partially 
agree with them; I can think of no good a priori reason to consider this particular quantity 
rather than some other. However, as long as (3 and W are both well-defined, then so is the 
value of e~@ w , and it is a perfectly legitimate exercise to investigate its properties. Any 
justification for embarking on such an investigation can only come a posteriori: if it leads 
us to an interesting and potentially useful result, then that is sufficient reason for studying 
it in the first place! 

My aim here is to address the arguments of CM, in particular their assertion that the 
derivation presented in Ref. is flawed. To do so, I will first present a detailed version of 
that derivation (Section H}, with the goal of establishing the nonequilbrium work relation 
as a mathematical identity, within the context of a Hamiltonian model for the evolution of 
the system of interest and its thermal environment. Following that, I will discuss several 
specific points raised by CM (Section ITTj) . 
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I. DERIVATION 



The nonequilibrium work relation (or "Jarzynski equality", in CM) can be stated as 
follows. Imagine a finite system that has been prepared in a state of equilibrium with 
a thermal environment at temperature T, and suppose that we subject this system to a 
thermodynamic process, by externally varying a work parameter of the system, A, from an 
initial value A to a final value B. In doing so we both drive the system out of equilibrium, 
and perform some amount of work, W, on it. The precise value of W depends of course on 
the specific motions of the microscopic degrees of freedom that constitute the system, and 
these motions are in turn influenced by the degrees of freedom of the environment. Therefore 
let us imagine that we carry out this process infinitely many times. During each of these 
repetitions, or realizations, of the process, we begin with the system and environment in a 
state of equilibrium, and we always vary the work parameter A in precisely the same manner 
from A to B. After each realization we note down the amount of work W performed on 
the system during that realization, and in the end we construct the distribution of work 
values, p{W), observed over this set of realizations of the process. The nonequilibrium work 
relation states that this distribution satisfies a strong constraint, namely: 

(e^ w ) = J dWp(W) e? w = e-? AF , (1) 

which remains valid regardless of how slowly or quickly we varied the work parameter during 
the process. Here, 1/(3 = ksT, and AF is the free energy difference between the equilibrium 
states associated with the initial and final values of the work parameter. To be precise, let 
Z\ denote the partition function (defined by Eq 1201 below) corresponding to the equilibrium 
state of the system of interest, when the work parameter is held fixed at the value A, and 
the system is in equilibrium with the environment. Then the free energy of that equilibrium 
state is given by the usual formula, F\ = — (3~ 1 \nZ\, and the quantity AF appearing in 
EqQ is defined to be 

AF = F B -F A = -/?- 1 ln|^. (2) 

A special case of Eq^ pertaining to the situation in which an external perturbation to 
the system is turned on and then off (hence AF = 0) was derived earlier by Bochkov and 
Kuzovlev. 

I will now give a detailed version of the derivation of the nonequilibrium work relation 
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found in Ref.[2j, and will begin by spelling out the assumptions behind this derivation. 

First, let us treat the system and its thermal environment as a set of classical degrees of 
freedom that are well isolated from the rest of the universe, and described by a Hamiltonian 



H(T; A) = H(x; A) + H E (y) + h- mt (x, y). 



(3) 



Here, x denotes a point in the phase space of the system of interest; y is a point in the 
(typically vastly larger) phase space of the thermal environment; T = (x, y) is a point in the 
combined phase space of system and environment; and A is an externally controlled work 
parameter. The terms on the right side of Eq|3] correspond to the bare Hamiltonian for the 
system of interest (if), the bare Hamiltonian for the environment (He), and the interaction 
energy between system and environment (/lint)- 

Let us now imagine that we prepare the system and environment in an initial state of ther- 
mal equilibrium at a temperature T, with the work parameter held fixed at an initial value 
A = A. To be specific, suppose that this preparation is accomplished by placing the combined 
system and environment in weak thermal contact with a much larger "super-environment" 
at temperature T, and then removing the super-environment after an appropriate equilibra- 
tion time. As a result of this preparation, the system and environment find themselves in a 
microstate T = (xo,yo) that can effectively be viewed as being sampled randomly from the 
canonical distribution in the full phase space: 



This is the classical partition function for the equilibrium state of the system and environ- 
ment, when the work parameter is held fixed at a value A. 

Having prepared the initial state of equilibrium and removed the super-environment, we 
allow the system and environment to evolve in time as we vary the work parameter from 
X = A at t = to X = B at t = r, according to some arbitrary but pre-determined schedule. 
The microscopic history of system and environment during this process is described by a 
trajectory {T t } evolving under Hamilton's equations in the full phase space. Here I use the 
notation {T t } to denote the entire trajectory, that is the microscopic history of the system 




(4) 



The normalization factor Ya is a particular case (A = A) of the quantity 




(5) 
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and environment from t = to t = r. By contrast, the notation Ft (without the braces) 
denotes simply the microstate of the system and environment at a specific time t. Similarly, 
{At} will refer to the schedule for varying the work parameter from A to B, and X t to the 
value of the work parameter at a particular time t. The schedule {X t } specifies how we act 
on the system during the process in question, whereas {T t } specifies how the system and 
environment respond, at the microscopic level, during a given realization of the process. 



Let us interpret H(x; A) as the internal energy of the system of interest. 17( The net 
change in this quantity during a single realization of the process is equal, identically, to 

H(x T ;B)-H(x Q ;A) = J dtX—{x t ,X t ) + J dtx—{x t ,X t ). (6) 

It is natural to interpret the first integral on the right side as the external work (or mechanical 
work in Ref.yJ) performed on the system, and the second term as the heat absorbed by the 



system (see for instance the discussions in Refs.P, 18, [3]). EqElcan then be viewed as a 
statement of the first law of thermodynamics, which asserts that the change in the internal 
energy of the system is due to two contributions: the work performed on the system, 

W = j o dtX—(x t ,Xt), (7) 

and the heat absorbed by the system, 

Q = j Q dtx—(x t ,X t ). (8) 

Eqs0 and |H1 define W and Q in terms of the microscopic history of the system alone, 
{xt}. In effect, while Nature integrates Hamilton's equations in the full phase space, we 
observers need to monitor only the degrees of freedom of the system of interest in order to 
deduce how much work was performed on the system, and how much heat was absorbed by 
it, during a given realization of the process. Note, however, that for a realization described 



by a trajectory {T t } in the full phase space, we have: [20] 



H(r T ;B)-n(T ;A) = £ dt jH(T t ; X t ) (9) 

dtx^(r t ;X t ) (io) 

OH 

dtX—(x t ;X t ), (11) 
o oX 

i.e. 

W = H(T T ;B)-H(T ;A). (12) 



This tells us that the work performed on the system of interest is equal to the net change 
in the Hamiltonian of the combined system and environment. 

The distinction between Eqs|7]andE3is an important one: the former defines the external 
work performed on the system, in terms of its microscopic evolution; the latter states that 
the quantity W thus defined is equal to the net change in the combined energy of the system 
and environment (under the assumption of Hamiltonian evolution for the combined system 
and environment). 

The preceding paragraphs have specified a model used to represent a system that is 
both in contact with a thermal environment, and subject to an externally controlled work 
parameter. Two strong but commonly made assumptions underlie this model: first, that 
quantum effects can be ignored; second, that the system and environment can be treated 
as being isolated from all other degrees of freedom. I will now show that Eq^ follows as a 
direct consequence of this model. 

So far the discussion has focused on a single realization of the thermodynamic process 
in question. Now imagine that we carry out the process very many - in principle, infinitely 
many - times. We always prepare the system and environment in equilibrium as described 
above, and we always use the same schedule {X t } to vary the work parameter from A 
to B. In other words, we act on the system in precisely the same way, over and over 
again. Nevertheless, the microscopic history of the system and environment, {T t }, will 
differ from one realization to the next, simply because the initial microstate To differs from 
one realization to the next (see EqE}. Over the course of each realization we observe 
the evolution of the system's degrees of freedom, {xt}, and from that empirical data we 
compute the value of W using EqJTJ Finally, from the set of work values collected over these 
realizations, we construct the average of e~^ w . Formally, 



where N denotes the number of realizations, and W n is the work performed on the system 
during the nth realization. As above, 1/(3 = kgT, where T is the initial temperature at 
which the system and environment are prepared, which is "the only known temperature 
available", as CM correctly point at the bottom of page 4 of Ref.|l|. 

The preceding paragraph describes how to construct the desired average from a series of 
measurements. Let us now analyze this quantity theoretically: for a given Hamiltonian 7i, 




(13) 
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initial temperature T, and schedule {Xt}, what value will we obtain for the average defined 
by Eqll3l.'' Note that this problem is fully specified: the probability distribution for the 
initial conditions T is given by EqHJ the subsequent evolution in the full phase space, {T t }, 
is determined by Hamilton's equations; and the quantities f3 and W are precisely defined. 
Thus the question we have just posed has a unique answer within the context of our model, 
and it remains only to do the math. 

To carry out the analysis, let us introduce a function W(Vq), which is the work performed 
on the system for a realization launched from initial conditions T in the full phase space. 
(The initial conditions T uniquely determine a trajectory {T t } in the full phase space, and 
from such a trajectory the value of W can be obtained through either EqJ7|or Eqll21) Then 
the average we wish to evaluate can be written as 



e -pw\ 



J dT p(T ) exp[-/W(r )], (14) 



with p(Tq) as given by Eq0] If we now combine Eq0] with Eqll2| then after a single 
cancellation we get 

(e-? w ) = J dT -L expH3W(r T (r ); B)], (15) 

where the notation stresses the fact that Hamiltonian evolution is deterministic: r r (ro) 
represents the final microstate in the full phase space, for the realization launched from 
the initial microstate IV Since there is a one-to-one correspondence between the initial 
conditions Tq and the final conditions T T , we can perform a change of variables in the above 
integral: 



J dT ■ ■ ■ = J dT T 



er • < 16 > 

where |<9r r /<9r | is the Jacobian associated with this change of variables. By Liouville's 



theorem, this Jacobian is equal to unity. 2l( We therefore have 



e -pw\ 



J dT T -L exp[-(3H(T T ; B)] = ^, (17) 

by EqUJ Since the ratio Yb/Ya depends only on the parameter values A = A and B, and 
on the temperature T, EqlTTI already establishes a strong result. Namely, even though the 
distribution of work values p(W) generally depends on the specific protocol for varying A 
from A to B, the quantity / dW p(W) e _/W does not! 

To this point, no approximations have been made in the analysis. Now, the right side 
of Eq|T7| is a ratio of partition functions of the combined system and environment (see 



EqISJ). We want to replace this by some expression pertaining to the system itself. In Ref. [3] 
this was accomplished by explicitly assuming the interaction e nerg y h- mt to be negligible 
in comparison with the other two terms in the Hamiltonian Ji. [22] In many situations of 
physical interest, however, h int is not negligibly small. In this case the evaluation of the ratio 
Yb/Ya requires a bit more effort. As a starting point, let us recall that when h- int is finite, 
the equilibrium distribution of the system of interest is given by the following modification 
of the familiar Boltzmann-Gibbs formula: 

p s (x; A) oc exp{-[3H*(x; A)], (18) 

where 

m , „ TT , Jdyexp[-p[H E (y) + h int (x,y)]] 
H*(x; A) = H(x; A) - (3 1 In J=- j - (19) 



I dyexp[-(3H E (y)} 

is a potential of mean force (PMF) associated with the phase space variables of the system of 
interest. [23] The notation ps indicates the probability distribution of the system of interest; 
this is obtained from the probability distribution in the full phase space, Eq01 by integrating 



24j For the equilibrium distribution given by 



over the environmental degrees of freedom. 
Eq J181 the partition function (normalization factor) is 

Z x = J dx exp{-(3H*(x; A)]. (20) 

With these definitions, we have 

Yx = Z x -J dy exp[-pH E (y)}, (21) 

which immediately implies 

Xl = ^. (22) 

Y A Z a K ' 

EqJ21lis not an approximation, but follows identically from Eqs|3E3 and l20l along with the 
definition of 7i, EqOl Combining these results with EqsEl and El we finally arrive at the 
nonequilibrium work relation: 

U-W\ = Yb = Zn = e _ pAF ^ (23) 

N ' i A ^A 

without resorting to a weak- coupling assumption. For physical situations in which hi nt hap- 
pens to be negligible, we recover the situation discussed in Ref. Q|. 
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Note that the quantity AF has been defined mathematically, in terms of a ratio of parti- 
tion functions fEqs l2l20"j) . In the Appendix, I briefly argue that it is reasonable to view this 
quantity physical free energy difference. 

It is important to stress that the derivation which has just been presented is exact: Eq 12*31 
is a mathematical equality, given the model specified above. The key feature of this model is 
that the system and environment are treated as an isolated collection of classical degrees of 
freedom evolving under Hamilton's equations. While this model represents the traditional 
approach of classical statistical mechanics, there is a subtlety associated with it, even if we 
agree to ignore quantum effects. This subtlety arises from the fact that it is in practice 
impossible to completely isolate a system and its immediate thermal environment; unavoid- 
able interactions with the rest of the universe introduce effectively random perturbations 
to the evolution of the trajectory {T t }. (Moreover these perturbations are correlated with 
{r t }.) No matter how weak these perturbations may be, their effect on a given trajectory 
becomes magnified exponentially with time if the evolution of {T t } is chaotic, as is typically 
the case for a realistic thermal environment. With this in mind, are we really justified in 
invoking Liouville's theorem in going from EqJ15lto EqU^li* This theorem reflects a delicate 
balance between sets of initial and final conditions of trajectories evolving under Hamilton's 
equations, a balance that might well be upset by the addition of even the tiniest amount 
of randomness into the equations of motion. Questions such as this make it all the more 
important that EqJ"Qbe tested in actual laboratory experiments. 

The above derivation relies explicitly on the assumption that the initial equilibrium state 
is represented by a canonical distribution in the full phase space (EqHJ). This assumption was 
justified by the somewhat artificial construct of a super-environment. However, the validity 
of Eq["[] might not depend as strongly on a literal interpretation of EqJfU as the derivation 
suggests. For instance, in the pulling experiment of Ref.jl^j], the microscopic system of 
interest - a single strand of RNA, two micron-size beads, and the DNA handles used to 
attached the RNA to the beads - is immersed in a macroscopic bath of water molecules. As 
long as that bath is prepared at a well defined temperature, one intuitively expects that the 
behavior of the biomolecule, immersed deep within the aqueous solution, will not depend in 
any significant way on whether the combined system and environment are prepared exactly 
in the canonical distribution given by Eq["f] Thus we would expect the same outcome (apart 
from extremely small corrections) if the initial conditions Tq were instead sampled from a 
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microcanonical distributi on, p oc 5(E — 7i); this is a variant of the usual expectation of 
equivalence of ensembles. [25| For a quantitative discussion of this issue, see Section II. B of 
Park and Schulten.ll! 



II. RESPONSE TO SPECIFIC POINTS 

A central claim of Ref . Q is that the heat exchange between the system and environment 
has "not been properly taken into account" (page 5 of Ref.jlJ.) While it is true that the 
derivation given above never makes explicit use of the quantity Q, this does not imply that 
Q is assumed to be zero, or that in some other way the heat has been mishandled in the 
analysis. It simply means that - within the context of the model - one can evaluate the 
average of exp(— f3W) without mentioning Q in the calculation. 

As an example of their claim that heat has not been treated properly, CM consider the 
situation in which the system of interest is out of equilibrium at the moment when the work 
parameter reaches its final value (page 5 of Ref. [l]). It is worthwhile to discuss this example 
in some detail. To begin, recall that AF should always be understood as the free energy 
difference between the two equilibrium states associated with the initial and final values of 
the work parameter, rather than as the free energy difference between the initial and final 
states of the system of interest. (Indeed, if the system is out of equilibrium at the end of 
the process, then its final free energy might not be well-defined.) To be absolutely precise, 
for any process during which the work parameter is varied from A to B, the quantity AF 
appearing in EqJTJis defined by (see Eqs J2l20|) 

= fdxex P {-(3H*(x;B)] 

Jdx exp[-(3H*(x;A)Y 1 1 

regardless of whether or not the system is in equilibrium at the end of the process. 

Now let us consider a two-stage schedule for varying the work parameter. During the 
first stage (0 < t < Ti), A is changed in some arbitrary way from A to B; during the second 
stage (n < t < T2), A is held fixed at the value B. Let us assume that the system is out of 
equilibrium at the time t%, but that the second ("relaxation") stage is sufficiently long for 
the system to relax to equilibrium. Now imagine two observers, one of whom monitors the 
behavior of the system only during the first stage, while the other monitors the behavior 
during both stages. Thus according to the first observer the process ends at time T\ (with 
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the system out of equilibrium), whereas the second observer contends that the process ends 
at a later time r 2 (with the system in equilibrium). As before, we imagine infinitely many 
repetitions of the process. 

For every realization, the two observers agree on the precise amount of work performed 
on the system during the process, even though they disagree as to when the process ends. 
This follows from Eq|7| since A is fixed for t > r%, there is no contribution to W from the 
relaxation stage. Therefore, when the two observers independently construct the average 
(exp(— (3W)) after many realizations of the process, they both arrive at the same number for 
the left side of Eq[0 Moreover, since both observers agree that the work parameter begins 
at A (at t = 0) and ends at B (at t = T\ or t = r 2 ), they also agree on the value of AF, 
as defined by Eql24l Hence when using their data to assess the validity of Eq^ the two 
observers will be comparing exactly the same numbers. In other words, whether or not we 
choose to include a relaxation stage - during which we hold A fixed at its final value, so as 
to let the system come to equilibrium - has no bearing whatsoever on the validity of EqO 
This is a simple consequence of the definitions of the quantities W and AF. 

Although no work is done on the system during the relaxation stage, there is typically a 
certain amount of heat exchanged between system and environment during this stage. Thus 
the observed values of Q generally do depend on whether the process is defined to end at 
time 7~i or time r 2 . But this in no way affects the validity of Eq^ since that prediction 
concerns work, as defined by Eq0 and not heat. 

Another point raised by CM concerns the factor (3 = l/fc^T, where T denotes the initial 
temperature of the system and environment. Once the system is driven away from equilib- 
rium, it might not have a well-defined temperature, and even if it did there is no guarantee 
that it would be equal to the initial temperature T. Therefore, in the expression e _/3M/ , a 
value that pertains to a system out of equilibrium (W) is divided by a temperature (T) 
that does not meaningfully represent the state of that system, except at t = 0. CM as- 
sert, with some justification, that such a pairing appears arbitrary and without foundation 
(see e.g. page 4 of Ref. [lj]). Note the nature of this criticism: CM do not claim that the 
quantity e~^ w is somehow inherently "illegal" or ill-defined, but rather that it is ad hoc. 
However, as discussed briefly in the introduction above, as long as (3 and W are well-defined, 
it is perfectly acceptable to inquire about the average of e~^ w over different realizations of 
the process. As the detailed calculation of Section H] reveals in the context of a particular 
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model - and as has been shown by a number of other derivations using significantly different 



models 
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ll\ - this average works out to be equal to e~ /3AF . It might be 
surprising that a construction as intuitively unnatural as e~^ w should lead to such a simple 
result, but this does not automatically invalidate the result. Indeed, the fact that this quan- 
tity seems unnatural might simply reflect our limited intuition regarding nonequilibrium 



processes. 



26] 



In Section 2 of their paper, CM consider two factors that play an important role in 
determining whether a process is reversible or irreversible. The first is the rate of heat 
transfer between the system and environment, c, which is related to the strength of the 
coupling between them; the second is the rate at which work is performed, w = XdH/dX. 
CM discuss several cases illustrating how the balance between c and w affects the reversibility 
or irreversibility of the process. Their discussion is physically motivated and certainly seems 
correct, but it does not bear on the validity of the derivation of the nonequilibrium work 
relation given in Ref.Q]. As shown in detail in Section |l] above, that derivation is based 
on very general considerations involving Hamilton's equations, Liuoville's theorem, and the 
use of the canonical ensemble to represent the initial equilibrium state of the system and 
environment. These considerations remain valid independently of the values of c and w. To 
put it another way: the analysis presented in Section |l] above depends neither on the rate at 
which the work parameter is varied, nor on the strength of the coupling between the system 
and its environment, hence it is as valid for irreversible processes as it is for reversible ones. 

In Section 5 of Ref . Q CM discuss the example of an isolated harmonic oscillator whose 
frequency is externally switched from an initial value uq to a final value u>i (where ui > ujq), 
over a switching time t s . In Ref. jj] I had shown that in the limit of infinitely slow switching, 
and assuming a canonical distribution of initial conditions for the oscillator, one can solve 
exactly for the distribution of work values: 



lim p{W) = ^-^ V (^^-\e(\\ ). (20) 



where #(■) is the unit step function. It is easy to verify that this distribution satisfies the 
nonequilibrium work relation, Eq^above. The derivation of EqJ2H]niakes use of an adiabatic 
invariant and therefore is valid only in the limit t s — > oo. This does not, however, imply that 
the nonequilibrium work relation fails for finite values of t s , only that for finite switching 
times it is not easy to obtain an exact expression for p(W). (Exactly solvable models are 
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hard to come by! See, however, Ref . [27 1 . ) Therefore in Ref. [3] the analytical calculation of 
p(W) for infinite switching times (Eq[2S]above) was supplemented by numerical experiments 
carried out for five different finite values of t s . The results of these experiments were in 
excellent agreement with EqQ as shown by the diamonds in Fig. 2 of Ref.jjJ. It is difficult 
to reconcile these results with CM's statement that the nonequilibrium work relation "is 
critically dependent on adiabatic invariance for finite t s , even for the harmonic oscillator" 
(page 14), particularly since in Ref.jjJ adiabatic invariance was only invoked in the limit 
t s — > 00. 

Of course a single example does not establish universal validity. General proofs of Eq^for 
Hamiltonian systems were given in Refs.j^j and jjj, and the harmonic oscillator was meant 
to serve only simple illustration. 

It is a pleasure to acknowledge useful discussions and correspondence with A. Adib, C. 
Bustamante, E.G.D. Cohen, G.E. Crooks, J. Jarzynski, J. Liphardt, and F. Ritort. This 
research was supported by the Department of Energy, under contract W-7405-ENG-36. 



Appendix 

Eql2*T1 can be written explicitly as 

J dTe-P*™ = J dxe-W**) . / dy C -/w-(»). (26) 

Thus the partition function for the combined system and environment factorizes nicely as 
the product of two partition functions, one for the system of interest (which includes all the 
effects of the interaction energy) and the other for the environment. If we take the natural 
logarithm of both sides and multiply by — we can rewrite the above result as 

— F\ + F°, (27) 

where T\ = —fl^ 1 In Ya can be viewed as the equilibrium free energy of the combined system 
and environment, and F% = — f3~ l In / dye~@ HE ^ as that of the bare environment. Note 
that is a macroscopic quantity, describing a macroscopic thermal environment, whereas 
the characteristic magnitude of Fx is determined by the size of the system of interest; e.g. 
Fa (and therefore AF) is microscopic for a single- molecule pulling experiment. 
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Given the definitions in Section HI it is easy to show that the quantity Fa = — (3 1 In Z\ 
satisfies dF x /d\ = (dH*/dX)^ = (dH/dX) e x q , or equivalently 

f B /dH\ eq 

Here (• • -)^ q = / dxps{x; A) • • • denotes an equilibrium average at a fixed value of the work 
parameter. (The derivation of EqEHJ not reproduced here, is just a few lines long, and 
essentially identical to the derivation of the well-known thermodynamic integration identity; 
see e.g. Ref.^.) 

Now suppose we carry out a reversible process: the system passes through a continuous 
sequence of equilibrium states as we slowly vary A from A to B. Thus, at any time t 
during this process, the system of interest is sampling its phase space according to the 
equilibrium distribution (ps) corresponding to the current value of work parameter, At. In 
this situation we can replace the value of dH/dX appearing in the definition of work (EqEJ), 
by its equilibrium average: 

-i-(f )>/>(!)>-< 

invoking Eql28l The arrow denotes that we are considering the special case of a reversible, 
quasi-static process. By these arguments, then, W = AF for any reversible process during 
which A is changed, quasi-statically, from A to B. This suggests that we are justified in 
interpreting AF, defined mathematically in terms of the modified Boltzmann distribution, 
Eq J181 as a physical equilibrium free energy difference. 
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